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Next generation aircraft with a large number of actuators will require advanced control 
allocation methods to compute the actuator commands needed to follow desired trajectories 
while respecting system constraints. Previously, algorithms were proposed to minimize the U 
or l 2 norms of the tracking error and of the actuator deflections. The paper discusses the 
alternative choice of the l x norm, or sup norm. Minimization of the control effort translates 
into the minimization of the ma xim um actuator deflection (min-max optimization). The 
paper shows how the problem can be solved effectively by converting it into a linear 
program and solving it using a simplex algorithm. Properties of the algorithm are also 
investigated through examples. In particular, the min-max criterion results in a type of load 
balancing, where the load is the desired command and the algorithm balances this load 
among various actuators. The solution using the l x norm also results in better robustness to 
failures and to lower sensitivity to nonlinearities in illustrative examples. 


I. Introduction 

Control allocation is the problem of distributing control effort among multiple, redundant actuators. In conventional 
flight control system design, the issue is resolved through the concept of ganging. Specifically, pseudo-effectors v 
are defined so that 

u = Gv (0.1) 


where v is the vector of pseudo-effectors, u is the vector of actuator commands, and G is a ganging matrix. For 
example, it is typical to define a single elevator command S e and a single aileron command 5 a so that 
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where 8 e] , 8 er are the left and right elevator commands, S a] , Sar are the left and right aileron commands, and 8 r is the 
rudder command. An elevator command 8 e produces symmetric left and right elevator commands, and an aileron 
command 8 a produces antisymmetric left and right aileron commands. This ganging typically results in a mostly 
decoupled response of the aircraft from the elevator, aileron, and rudder commands to the pitch, roll, and yaw 
responses. 

In control design for future vehicles, reasons to look for alternatives to ganging include cases where: 

1) the vehicle has a large number of actuators, making it less intuitive how the ganging matrix should be 
defined; 
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2) the vehicle has unconventional control effectors, with a significant degree of nonlinearity and of interaction 

between the effectors, again making it difficult to develop an intuitive solution; 

3) the effectiveness of the actuators is limited, making it important to optimize their use within their position 

and rate limits; 

4) the control system is designed to be reconfigurable or adaptive, so that control allocation must he computed 

in real-time. 

Previous work in control allocation includes the seminal paper of Durham 11 , which introduced the concept of 
direct allocation. Although the concept was mathematically formulated, algorithms to solve the problem numerically 
in the general case were only later developed 12 . In the meantime, Buffington 2 proposed an alternative formulation 
minimizing the norm of the error between desired and achieved commands. Using the h norm, he showed how the 
problem could be converted to a linear program and solved exactly, using standard linear programming software. 
Ikeda and Hood 18 similarly reported the application of U optimization, although with fewer details. Nevertheless, it 
became clear that the solution of optimal control allocation problems was feasible in real-time. In Ref. 22, it was 
shown that the direct allocation problem could also be solved using linear programming, and that a considerably 
smaller linear program could be obtained for the U optimization problem, compared to Ref. 3. Timing data showed 
that solutions of the problem could comfortably be performed in real-time, even for large numbers of actuators, and 
that the optimal solution improved performance significantly over simpler, approximate methods. 

Solutions of the optimal control allocation problems using the l 2 norm were also proposed, with an early solution 
provided through the fixed-point method of Ref. 4. The fixed-point algorithm was extremely simple, and many of 
the computations needed to be performed only once, before iterations started. Remarkably, the algorithm also 
provided an exact solution to the optimization problem and was guaranteed to converge. Numerical tests, however, 
showed convergence of the algorithm could be very slow and strongly depended on the problem (the number of 
iterations required could vary by orders of magnitude, depending on the desired command). An elegant alternative to 
this algorithm was proposed by Harkegard, using the theory of active sets 17 . The algorithm was very similar to the 
simplex algorithm used for l, optimization, and had the same advantage of completing in finite time and with a small 
number of iterations. 

Interior-point methods were also studied to solve large control allocation problems, both for the l, norm 21 and for 
the h norm 22 . The computational requirements of these methods scaled better with the number of actuators, but the 
number of actuators had to be quite large (>1 5) before the advantages become apparent. 

Among recent work, one may note several papers considering the application of control allocation to hypersonic 
vehicles, including some flight tests 9 ' 10, 23 , the problems posed by nonlinear actuator effectiveness (e.g., Ref. 8), 
modifications to account for the dynamic response of the actuators (e.g., Ref. 20), and the combination of control 
allocation with adaptation (e.g., Ref. 25, 26). New aircraft, for example those with blended wing body 
configurations 7 have been identified as presenting control allocation challenges due to novel actuators, distributed 
actuators with low control authority, interactions between control effectors, and interactions between propulsion and 
control surfaces. Interestingly, control allocation is also emerging in other applications, including land and marine 
vehicles 13, 15,24 . 


II. Optimization formulations of control allocation 


A. Control allocation in model reference control 

We introduce control allocation in the context of model reference control (a form of dynamic inversion). However, 
solutions may be used in a variety of control design methods. To state the problem mathematically, we consider the 
state-space model 


x A = A a x a +Bu + d 
yA = Cx a 


(0.3) 


where x A £R",(/6 R n , u E R p , y A E R q . For the control of aircraft, the states are given by the vector x A and may 
include the angle of attack, the pitch rate, the angle of sideslip, the roll rate, and the yaw rate (n=S). The output 
vector^ may contain the pitch rate, the roll rate, and the yaw rate (q=3). The control input vector u consists of the 
commanded actuator positions. In a conventional aircraft, these commands are the deflections of the two elevators, 
the two ailerons, and the rudder (p=5). The disturbance vector d represents the forces and moments that the control 
surfaces must cancel in order to trim the aircraft (i.e., to create an equilibrium of the dynamical system). 

For the purpose of example, consider a simple model reference control law. The method relies on a reference 
model that represents the desired dynamics of the closed-loop system 
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■>V - + Bm T M (0.4) 

where r M is a reference input vector (the pilot commands) and y M represents the desired output of the system. Since 
the derivative of y is given by 

y A = CA a x a + CBu + Cd (0.5) 

the objective may be achieved by setting 

CBu = -CA a x a -Cd + A M y A + B M r M = a d (0.6) 

where a d represents the desired vector to be matched by CBu. If y is a vector composed of the rotational rates (as is 
typically the case), a d represents desired rotational accelerations. 

Obtaining u from a d requires that one solve a system of linear equations with more unknowns than equations. 
Solving such a system is easy, but the difficulty in control allocation is that the vector u is constrained. The limits 
generally have the form 

“min,; £ K, s “max,/ f°r *' = 1 ,-,P (0.7) 

or, w min <.u<. « max , in vector form. There may be additional constraints due to the maximum rate of deflection of 
the actuators. We refer to the problem of finding a vector u that is the “best” possible solution of (0.6) within the 
constraints (0.7) as the control allocation problem. 

Given the constraints, the control allocation problem may be such that: 

• many solutions exist, 

• only one solution exists, 

• no exact solution exists. 

One is naturally drawn to finding solutions that minimize the error CBu-a d . Indeed, providing all the control 
authority available may make the difference between a maneuver being achievable or not, and between an unusual 
condition being recoverable from or not. However, the question also arises as to which solution is the most desirable 
when many solutions exist. Therefore, control allocation typically consists both in error minimization and control 
optimization. 

B. Formulations of control allocation 

The fundamental control allocation problem can be formulated as the following error minimization objective. 

Error minimization: given a matrix CB, find a vector u such that 

J = \CBu - a d \ (0.8) 


is minimized, subject to w min sms w max . 

The problem is solved exactly if 7=0. However, regardless of whether an exact solution exists, the following control 
minimization problem may be considered as well. 

Control minimization: given a matrix CB, a vector u p , and a vector u x such that w mjn s u x s « max , find a vector u 
such that 


J = u -u p 


(0.9) 


is minimized, subject to 

(CB)u = (C5)n, 


( 0 . 10 ) 


and u mm sms M max . 


The control minimization problem is a secondary optimization objective to be satisfied if the solution of the primary 
objective, given by u x , is not unique. The vector u p represents some preferred position of the actuators (e.g., zero 


deflections). After a solution yielding minimum error is obtained, the solution with minimum deviation from the 
preferred position is picked among all equivalent solutions. For both problems, weighting of the elements of the 
vectors may be inserted in the norms, either to prioritize the axes or to prioritize the actuators. 

The norm used in the optimization criteria is a design choice that has more consequences than might be 
expected. The h norm of a vector x is the sum of the absolute values of the elements of the vector 
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( 0 . 11 ) 


14 = 2 1** I 

i - 1 

while the h norm is the usual Euclidean norm 


( 0 . 12 ) 

Algorithms have been proposed for both norms and the results of the optimization problems are sometimes quite 
different. 

A possible implementation of optimization for control allocation consists in the sequential minimization of the 
error vector and of the control vector. Specifically, the error is minimized first, and then the control vector is 
minimized among all equivalent solutions. In Ref. 3, the control minimization problem was solved only when the 
solution of the primary error minimization problem was J= 0. However, it should be noted that, unless the matrix CB 
satisfies specific conditions (any q*q submatrix of CB must be nonsingular), the solution is not necessarily unique, 
even if the desired vector a d is not feasible. Given this fact, mixed optimization makes sense, and has several 
advantages over sequential optimization. 

Mixed optimization: Given a matrix CB and a vector u p , find a vector u such that 

J=\CBu-a d \+s\u-u p \ (0.13) 

is minimized, subject to u min s u s u max . 

The mixed optimization problem combines the error and control minimization problems into a single problem 
through the use of a small parameter £ . If the parameter £ is small, priority is given to error minimization over 
control minimization, as is normally desired. Often, the combined problem may be solved faster, and with better 
numerical properties, than when the error and control minimization problems are solved sequentially 1 . 



C. Optimization using the l, norm 

In this section, we review how the mixed lj optimization problem can be converted to a linear program of small 
size, following the presentation of Ref. 1. Further derivations later in the paper will build on this background. 

A standard linear programming problem consists of finding a vector x such that 

J = c T x (0.14) 

is minimized, subject to 

Osxsh, and Ax = b (0.15) 

In (0.15), vector inequalities are to be interpreted element-by-element. Alternative formulations exist, replacing 
Osrs/i by x a 0 , and Ax = b by Ax a b . However, these differences are not significant and the present form is 
preferable for the control allocation problem. 

For the conversion of the mixed optimization problem, define the function ,s(x) 

s(x) = x if x > 0 

' A . (0.16) 

= 0 otherwise 


This function is to be interpreted element-by-element in the vector case. We assume that the preferred vector 
satisfies u min u max . This condition may be eliminated without much difficulty, once the technique is 

understood. Define 

u + =s(u-u p ), u~=-s(u p -u ) (0.17) 


so that 

m = m + -u~ +u p 

0 * u+ s «max -«p> 0 S u ~ £ Up -Umin 

Similarly, define 

e = CBu -a d , e + = s(e), e~ = -s(-e) 


(0.18) 

(0.19) 
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so that 


Q e e , 0 £ Q £ ^max 5 0 £ ^ ^ ^rr 


( 0 . 20 ) 


where e max is some upper bound on the achievable error, e.g., e max = CBu p - a d | . 

With these definitions, the optimization problem involves a system of n linear equations 
e + - e~ - CBu + + CBu~ = CBu p - a d 

and the cost criterion 

• /= 5X + 5X +£ E m ' + +£ 2 m ' 7 

i-l M i=l /=1 

Therefore, defining the vector x T = (e + e~ u + w~), the linear programming problem is specified by 
A = (l -I -CB CB ), b = CBu p - a d 


( 0 . 21 ) 


( 0 . 22 ) 


-(1 


1 e 


0 


(0.23) 


hT = (<W ®max “max " U p u p ~ “min ) 


Note that the vector c implements an equal weighting of the elements within the vector a d and within u. However, 
the elements of the vector c can be changed to account for various objectives. For example, Ref. 3 showed how 
various choices could be made to minimize drag, wing loading, radar signature, or the use of thrust vectoring. 
Interestingly, different weights can even be applied for positive and negative values. 

Note that the A matrix of the linear programming problem has as many rows as the CB matrix. For the standard 
case with a 3-dimensional vector a d , the number of rows is only 3. This size is very small in linear programming, so 
the problem can be solved in a few iterations using, for example, the simplex algorithm. The algorithm is guaranteed 
to find an optimal solution in a finite period of time, it is easy to code, and it works well in practice. Speed of 
algorithm execution can be minimized by taking advantage of particular aspects of the control allocation problem. 
Because the number of columns in the A matrix is typically much greater than the number of rows, the problem is 
well suited for the so-called revised simplex method. The number of computations in this method depends only 
moderately on the number of columns. Because most variables of the vector x naturally have both upper and lower 
bounds, it is also advantageous to implement a simplex algorithm with both bounds, as opposed to the more 
common method with lower bounds which requires a large number of so-called slack variables. 

Speed of execution may also be improved significantly by initializing the simplex algorithm with a so-called 
basic feasible solution. Without this special feature, an initialization phase has to be added, requiring the use of 
another implementation of the simplex algorithm. A basic feasible solution is a vector x that solves Ax=b and is such 
that all elements of x are at their limits except q elements, where q is the number of rows of A. The mixed 
optimization algorithm can be initialized with u=u p as a feasible solution, so that 


u + =0, e + = s(CBu p -a d ) 
u =0, e~ = -s(-CBu p +a d ) 


(0.24) 


In general, q elements of e and e will be equal to zero, leaving only q elements (or fewer) different from zero. 
These elements are the basic variables of the initial feasible solution. Other usefiil techniques for the implementation 
of the revised simplex algorithm include the Sherman-Morrison-Woodbury formula to reduce the size of the 
matrices to be inverted and anticycling procedures to avoid infinite loops 1 . 


D. Implementation of control allocation solutions 

The implementation of optimal control allocation methods in real-time has become feasible. It is likely that obstacles 
will lie solely with validation and certification issues 1 . Generally, active sets and simplex methods require a finite 
number of steps for convergence, but the theoretical maximum is significantly greater than the typical number 
required. Further, the theoretical maximum assumes perfect computations. In the near term, the most viable 
implementation of optimal control allocation may be in the form of table look-up 14 . Control allocation may also be 
implemented as an add-on module that is engaged when failures are detected. 
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An issue for implementation of control allocation algorithms is the rate limits of the actuators. These can be 
accommodated by reducing the limits u m m, u max by the amount imposed by the rate limits. The preference vector 
may also be set to be the current actuator position. In this paper, we do not consider rate limits, focusing instead on 
the properties of control allocation solutions in general terms. 


III. Control allocation with load balancing 


A. Properties of /, optimization problems and load balancing 

Linear programming theory implies certain properties of the solution of the mixed /; optimization problem. 
Specifically, if the matrix CB has 3 rows, all the elements of the optimal vector x except 3 will be either at then- 
upper limit or at their lower limit. In terms of the control vector, this property implies that all but three control 
variables will be either at the upper limit, at the lower limit, or at the preferred position. If the vector a d cannot be 
achieved in any direction, all the control variables will be at one of the limits or at the preferred positions. The 
desirability of this property may be debated: on the one hand, it makes sense if the algorithm does not use ineffective 
surfaces. On the other hand, it is desirable to see all surfaces move together to achieve the desired moment. A more 
balanced distribution of the required effort to the control surfaces reduces the chances of encountering the control 
surface rate limits. 

In this section, we consider control minimization using the l x norm 

[lull = max \u t \ (0.25) 

00 i 

The norm of a vector is the maximum of the absolute values of the elements of the vector. It is also called the sup 
norm. For control optimization, use of the l„ norm in 

J = \\u-u II (0.26) 

II * I loo 

leads to an optimization criterion referred to as a min-max criterion, since the objective is to minimize the maximum 
value (in absolute terms) of the elements of the vector. This criterion has been used in a variety of networking 
control problems, including communication networks 19 and computer networks 16 . Typically, this criterion arises 
when attempting to balance the loads among multiple resources, such as processors or communication nodes. The 
or min-max criterion has also been used in control design 5 , although not as frequently as the l 2 criterion. 

In control allocation, use of the L, norm implies that one attempts to minimize the deflection of the actuators in 
the min-max sense. It does not matter how many actuators move: the maximum deflection should just be as small as 
possible. The solution that is obtained reflects this choice, by providing a more balanced distribution of the 
deflections than with the l t norm. Interestingly, the control allocation problems using the /^norm can be converted 
to linear programs that are similar to the l t linear programs, and solved using the same algorithms. 


B. Mixed //-/«, optimization 

We first consider the optimization of the criterion 

J=||CBM-a rf |[ 1 +e||w- (0.27) 


In other words, the // norm is used for the error minimization and the l x norm is used for control minimization, with 
both criteria mixed in a single, mixed optimization criterion. A small modification of the approach used for mixed h 
optimization yields the desired linear program. 

Introduce an additional variable u, which is intended to become the l m norm of u-u p . Next, vectors of slack 
variables 8« + and bu are introduced such that 


du + =u -u + 
8u~ = u -u~ 


(0.28) 


Using the same notation as for the /; optimization, a linear program can be defined with the cost criterion 

J = '2 e ? + ^ e 7 + 

;=1 i = 1 

and the optimization vector 


(0.29) 
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x T = (e + e u + 
The linear program to be solved is 


Su + Su 
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lx/; ^ 1 X p 


(0.30) 


(0.31) 


k T -(e, 


“mn-Up 


Up w m:n 


* \ 

Up ^min ^max I 


where I a «b is the identity matrix of dimension axb, 0 a «b is a matrix of dimension axb filled with zeros, l a «b is a matrix 
of dimension axb filled with ones, and 


"max =maX^« max -Mj 


\\ U p- U m in 


(0.32) 


Since the mixed // -/«, control allocation problem can be converted to a linear program, standard algorithms can 
be applied to solve it efficiently. Initialization with a basic feasible solution can be performed as for the // 
optimization, by adding to the original basic variables the new variables 8 k + and Sh\ The major drawback is that the 
number of rows has grown considerably in the process. From q rows (typically 3), the number has grown to q+2p 
(where p is the number of actuators). Nevertheless, such problems can still be solved very quickly on standard 
computing hardware. 


C. Mixed /oo optimization 

The mixed /«, control allocation problem is defined by the criterion 

y=|CSM-a rf | oj +e||M-M i ,|| (0.33) 


This problem can be converted to a linear program using similar techniques as the mixed // -l m control allocation 
problem. Specifically, define an additional variable e, which is intended to become the /« norm of e. Next, 
introduce vectors of slack variables Se + and be such that 


8e + =e‘-e + 


8e 


= e 


-e 


(0.34) 


A linear program can be defined with the cost criterion 

J = e + eu (0.35) 

and the optimization vector 

x T = (e + e~ u + u~ 8u + du~ u 8e* bu~ e ) (0.36) 

The linear program to be solved is 
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^ V^max e max ^max ^min ^max ^ p ^p ^min ^max ^max ^max ^max J 

An issue with this new formulation is that the computations grow even more, due to the further increase in the 
number of rows. It is possible that a smarter implementation will produce a more efficient algorithm. For example, 
Ref. 6 has an algorithm that is said to be more efficient than the simplex algorithm. However, constraints are not 
included. 

Our evaluation of the mixed l x criterion in specific examples has not shown significant differences with the 
mixed l r L optimization criterion. One reason is that any difference can only be noticed for non-feasible acceleration 
vectors (an achievable vector yields a zero error no matter what norm is used). Another reason is that the error 
vector has a small dimension compared to the control vector, and differences between the norms only become 
apparent when the dimensions of the vectors are large. 


IV. Numerical Results 


A. Low dimension example 

We first consider the aircraft model used by Durham 11 , which comes from NASA Dryden’s “Controls Design 
Challenge”, for a flight condition at Mach 0.5, 10,000 ft altitude. The CB matrix and actuator limits are given by 


/ 7.35 7.55 -1.35' 

[o.856 5.13 -13.7 y 


(0.38) 


Umax =[20,20,30]' 

^min ^max 


(0.39) 


The original CB matrix was multiplied here by 10 4 for numerical reasons, and the outputs associated with CB are 
moments, instead of accelerations, but these differences are insignificant for this example. The rows of the CB 
matrix are associated with roll and yaw moments, and the commands are the ailerons, the differential horizontal tail, 
and the rudder. 


Fig. 1 shows the results of three control allocation algorithms for a roll command. The x-axis gives the roll 
command, as a percentage of the maximum achievable pure roll command. The top plot shows the aileron 
command, and the bottom plot shows the differential horizontal tail command. The plots show the results of mixed 
l, optimization using the algorithm of Ref. 1 (solid), of mixed l 2 optimization using the algorithm of Ref. 17 
(dashed), and of mixed l,-L optimization using the algorithm of Ref. 1 modified as indicated in this paper (dotted). 
For mixed l 2 optimization, the criterion that is optimized is 

J = \\CBu - a d f 2 + s 2 \\uf 2 (0.40) 


subject to u mi „ sms w max . Note that the norms are squared, as well as the control weighting parameter e. The 

parameter e was set to 10' 3 for mixed l, optimization and mixed l r l„ optimization, and 10' 6 for mixed l 2 
optimization. A smaller parameter was set in the l 2 optimization to avoid solutions with acceleration errors. 
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Figure 1. Aileron and differential horizontal tail commands as functions of roll acceleration 
(solid: l t optimization, dashed: l 2 optimization, dotted: /,-/«, optimization). 

The top plot in Fig. 1 shows that the mixed l, optimization uses the most effective surface (the ailerons) to 
achieve the desired result for commands up to 50%. Once the limit is reached, the second most effective surface (the 
differential tail) is used. Note that the differential horizontal tail is actually slightly more effective at producing a roll 
moment. However, the significant yaw moment it produces makes it less effective than the aileron in the sense 
specified by the optimization criterion. The mixed l 2 optimization and mixed l r L optimization produce similar 
results, using both actuators from the beginning of the plot to the end. However, the mixed /;-/«, optimization 
produces equal deflections of the two actuators. In this case, optimal load balancing results in load equalization (this 
is not true in general for control allocation). 

As pointed out in Ref. 14, the fact that mixed h optimization only uses the most effective actuator may cause 
problems if the main actuator for a given axis fails. Fig. 2 shows the acceleration that is achieved in the case of an 
aileron failure (top) and differential horizontal tail failure (bottom). Since only the ailerons are used in the case of a 
small command, the aileron failure results in a loss of response when using the mixed // optimization. In fact, a 
small negative response is observed, due to negative moment produced by the rudder (which is applied to 
compensate for the yaw moment that would arise if the ailerons were effective). For larger commands, the system 
responds when the differential horizontal tail kicks in. 

Conversely, a failure of the differential tail does not degrade the performance of the system with mixed /; 
optimization for small roll commands, because the actuator is not used. Degradation is found when the commands 
become large. The mixed l 2 optimization and the mixed Ij-h, optimization degrade more gracefully, reducing the 
response for either failure, but maintaining a monotonically growing response throughout the range. The mixed /;-/«, 
optimization actually maintains a linear response, although with reduced magnitude, throughout the range of 
operation and for both failures.. 
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Figure 2. Roll acceleration obtained after failure of aileron (top) and differential tail 
(bottom) (solid: l, optimization, dashed: l 2 optimization, dotted: l,-l K optimization). 


B. Tailless aircraft example 

The second example is based on Lockheed's ICE ( innovative control effectors) tailless aircraft model found in 
Ref. 2. It has 11 actuators: left elevon, right elevon, pitch flaps, left all-moving tip, right all-moving tip, pitch thrust 
vectoring, yaw thrust vectoring, left spoiler slots, right spoiler slots, left outboard leading-edge flaps, and right 
outboard leading-edge flaps. The CB matrix corresponding to an output vector composed of pitch rate, roll rate, and 
yaw rate is given in Ref. 2 as 


/ -2.51 14 

-2.5115 

-1.9042 

-0.9494 -0.9494 


3.7830 

-3.7830 

0 

1.8255 

-1.8255 


v 0.0453 

-0.0453 

0 

-0.2081 

0.2081 


-1.1329 

0 

1.5046 

1.5046 -0.0003 

-0.0004 \ 

0 

0.0790 

-2.0956 

2.0957 -0.3067 

0.3067 

0 

-0.8038 

-0.0283 

0.0283 

0.0937 

-0.0937 / 


(0.41) 


for a flight condition at Mach 0.4 and 15,000 ft altitude. The position limits are given in Ref. 3 as 
"max = (30 30 30 60 60 10 10 10 10 40 40)' 


(0.42) 


M min “ (“30 -30 -30 0 0 -10 -10 0 0 0 0)' 


(0.43) 


Note that the limits of the spoiler slot deflectors were lowered from 60 degrees to 10 degrees in Ref. 3 to reduce 
nonlinear interactions between the spoiler slot deflectors and the elevons. The same limits were used here, although 
the nonlinear effects were not part of the evaluation. 

Fig. 3 shows the //, l 2 , and lx norms of the control vector, as functions of the percentage of maximum yaw 
acceleration. As may be expected, the l] norm is minimized for the mixed // optimization, the h norm is minimized 
for the mixed l 2 optimization, and the lx norm is minimized for the mixed l r lx optimization. This is generally true, 
but not always, due to the mixed nature of the optimization criteria. In Fig. 3, one finds that the growth of the lx 
norm of the control vector is considerably delayed with the mixed h- lx optimization compared to the other 
methods, even though the h and h norms are quite comparable. The extra room that is provided away from the 
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limits may be useful, for example, to add additional excitation for the purpose of real-time parameter identification 
(e.g., the null space injection discussed in Ref. 3). 

Fig. 4 shows how the reduction in the magnitude of the control vector with the mixed l r lx optimization may 
yield better performance, by avoiding end-of-range nonlinearities in the actuator effectiveness. Specifically, the 
figures assume that the actual acceleration vector produced is of the form 

a = Cb(u- au 3 ) (0.44) 

In other words, a small cubic nonlinearity reduces the actuators effectiveness uniformly as the deflections increase. 
The function u 3 is to be interpreted element by element. The constant a was set so that the effectiveness of a given 
element was reduced by 5% when the magnitude of the element of u was 60 degrees. This is a fairly small amount of 
nonlinearity, especially since most surfaces are limited well before 60 degrees. 

The nonlinearity only yields to a reduction of yaw acceleration of the order of 3% at the end of the range. 
However, note that a significant cross-coupling error occurs in the roll axis. In all plots, one finds that the mixed l r 
loo optimization criterion yields the lowest error, by virtue of the smaller actuator commands that result from its 
solution. 



Figure 3. l h l 2 , and /«> norms of the control vector, as functions of the percentage of 
maximum yaw acceleration (solid: h optimization, dashed: l 2 optimization, dotted: l r lx 
optimization). 
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Figure 4. Yaw and roll errors due to nonlinear actuator response as functions of the 
percentage of maximum yaw acceleration (solid: lj optimization, dashed: l 2 optimization, 
dotted: /,-/«, optimization). 


V. Conclusions 

The paper first reviewed how the mixed l, optimization criterion could be converted into a linear program and 
solved using a simplex algorithm. Then, it was shown that the formulation could be modified to solve mixed l r lx (// 
for error, /«, for control) and mixed l x (/« for error and control) optimization problems. It was argued that the l x 
norm leads to better “load balancing”, as defined in networks. A numerical example indeed showed a better balance 
in the use of the actuators. Advantages of such a feature were shown to include a greater resilience to actuator 
failures and to nonlinear effectiveness for large actuator deflections. The goal of the paper was not to prove that one 
method is better than another: rather, it was to increase the number of choices available to the engineer, as well as 
the understanding of how the choices relate to important properties of the solutions. 
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